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Within the framework of macroscopic quantum electrodynamics, general expressions for the 
Casimir force acting on linearly and causally responding magnetodielectric bodies that can be em- 
bedded in another linear and causal magnetodielectric medium are derived. Consistency with micro- 
scopic harmonic-oscillator models of the matter is shown. The theory is applied to planar structures 
and proper generalizations of Casimir's and Lifshitz-type formulas are given. 
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I. INTRODUCTION 

It is well known that the introduction of momentum 
and energy of the macroscopic electromagnetic field re- 
quires careful consideration, even for linear media. In 
fact, insertion of constitutive relations into Gauss's and 
Ampere's laws prevents one, in general, from deriving 
local balance equations of a similar type as in the micro- 
scopic theory. However, with quite restrictive (and most 
questionable) assumptions about the material under con- 
sideration, these difficulties can be formally overcome. 
Therefore, textbooks typically resort to approximate for- 
mulas that are based on assumptions such as quasi- 
monochromatic fields and lossless media [l|, |2J . Although 
the limitations inherent in such theories are rather obvi- 
ous, they are nevertheless applied beyond their range of 
validity. 

A typical example is the Casimir effect, which is closely 
related to the changes in the vacuum electromagnetic- 
field energy and/or momentum flow (stress) induced by 
the presence of inhomogeneous matter. With regard to 
the calculation of the Casimir force acting on macroscopic 
bodies that are embedded in a medium, the question of 
what are the correct expressions for these quantities be- 
comes crucial. Frequently, expressions that seem reason- 
able at first glance — such as Minkowski's stress tensor — 
have been taken for granted without justification. As we 
shall see, this has led to incorrect extensions of the well- 
known Lifshitz formula for the Casimir force between two 
dielectric half-spaces separated by vacuum to the case 
where the interspace is not empty but also filled with 
material 0, 0| ( see also the textbooks 0, 0, 0] and 
references therein). 

In this paper, we reconsider, within the framework 
of macroscopic quantum electrodynamics, the problem 
of the calculation of Casimir forces, by regarding the 
Lorentz force density as the fundamental quantity. The 
Lorentz force acting on some (macroscopic) spatial re- 
gion containing (electrically neutral) matter is of course 
the corresponding volume integral of the Lorentz force 
density, where the relevant charge and current densities 
may be thought of as being expressed in terms of the po- 
larization and the magnetization of the matter. As a con- 
sequence, the Casimir force on a macroscopic body and, 



equivalently, the stress on its surface can be expressed — 
in close analogy with microscopic electrodynamics — in 
terms of the electric and induction fields, irrespective of 
any specific constitutive relations. In particular, if the 
body linearly responds to the electric and induction fields 
and the medium the body is embedded in is also a lin- 
ear one, then the Casimir force can be expressed solely 
in terms of the classical (retarded) Green tensor, which 
in turn is determined by the response functions of the 
magnetodielectric matter under consideration. 

We show that the Casimir force formula found in this 
way is consistent with microscopic theories based on 
harmonic-oscillator models of dispersing and absorbing 
dielectric matter. The formula is valid under very gen- 
eral conditions, and enables one to study in a consistent 
way the Casimir force on linearly responding, dispersing, 
and absorbing magnetodielectric bodies that are not nec- 
essarily placed in vacuum but may also be surrounded 
by a dispersing and absorbing linear magnetodielectric 
medium. Since magnetodielectric matter is included in 
the theory, it is possible to consider also left-handed ma- 
terials 0| . To illustrate the theory, we apply it to a pla- 
nar geometry and derive a proper extension of Lifshitz- 
type formulas, with emphasis also on the extension of 
Casimir's original formula. 

The paper is organized as follows. In Sec. [HJ the 
stress tensor associated with the (macroscopic) Lorentz 
force is introduced. The Casimir force is calculated in 
Sec. IIIII and Sec. lIVI makes contact with the microscopic 
harmonic-oscillator model. The application of the theory 
to planar structures is given in Sec. [V] and a summary 
and some concluding remarks are given in Sec. IVII 



II. LORENTZ FORCE AND STRESS TENSOR 

Let us begin with the classical Maxwell equations for 
the electric and induction fields E and B in the presence 
of matter, 
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In this paper, "dot products" are written without the 
dot, and dyadic products are denoted by (g>. In Eqs. © 
and p and j cover all charges and currents of the 
system under consideration. Within the framework of a 
macroscopic description, the "internal" charges and cur- 
rents associated with the particles that form some neu- 
tral material system are commonly described in terms of 
polarization and magnetization fields P and M, respec- 
tively. The remaining "external" charges and currents — 
if any — are kept explicitly, i.e., 



P = Pint + Pcxt, 
J Jint ~t~ Jcxt ) 



where 
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As long as constitutive equations (relating the polariza- 
tion and magnetization fields to the electric and induc- 
tion fields) are not introduced, Eqs. which may 
also be interpreted microscopically, are generally valid. 
Clearly, the defining equations Q and JSJ of P and M, 
respectively, can be satisfied for any choice of (conserved) 
"internal" sources, the corresponding integrability condi- 
tion being just 



dpu 



dt 
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Note that the "internal" sources typically comprise 
bound charges and the associated currents — a concept 
that is commonly used together with spatial averaging in 
macroscopic electrodynamics [Io|. 
As known, the Lorentz force density 

f = pE+jxB (10) 

can be rewritten with the help of Eqs. ©-(0} as 



f = VT-eo^(ExB), 



(11) 



where the stress tensor 



T = e E <g) E + ii^B OB- ±(e E 2 + ^B 2 )l (12) 

has been introduced (1 , unit tensor). Clearly, Eqs. i|10|) 
and 1)11(1 are universally valid, regardless of whether the 
charge and current densities have been decomposed ac- 
cording to Eqs. ©-© or not. Note that essentially 
the same position has been recently taken up |ll| in 
the (re)analysis of measurements of the electromagnetic 
force that acts on dielectric 0, ITM Jbfl. ITEl IT^ | or mag- 
netodielectric 01 ( see & l so Refs. jl8Hl9f) disks exposed 



to crossed electric and magnetic fields. (For a different 
perspective, see also Ref. |20| ). 

The integral of the Lorentz force density f over some 
space region (volume V) gives of course the total electro- 
magnetic force F acting on the matter inside it, 



F= / d 3 rf. 
Iv 



Using Eq. (fTTJl . we have 
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which is obviously also true if the space region is occupied 
by a macroscopic body, with the charges and currents be- 
ing "internal" ones described by polarization and mag- 
netization fields. In particular, if the volume integral on 
the right-hand side of this equation does not depend on 
time, then the total force reduces to the surface integral 



dF, 



where 



dF = daT = Tda 



(15) 



(16) 



may be regarded as the infinitesimal force clement act- 
ing on an infinitesimal surface element da. Note that a 
constant term in the stress tensor does not contribute to 
the integral in Eq. I|15|l and can therefore be omitted. In 
the calculation of the Casimir force in Sec. II I II it will be 
necessary to make use of this fact. 

Expressing in Eq. I|16f) the stress tensor T in terms of 
Minkowski's stress tensor T*- M ^ (which agrees with Abra- 
ham's stress tensor 0), 

T (M) = D<g)E + H<g)B — i (DE + HB) 1 

= T + P®E-M®B-|(PE-MB)1, (17) 

one finds that 

dF = daT (M) -da[P®E-M®B - |(PE - MB)] , 
from which it is seen that in general 



dF ^ daT 



(M) 



(19) 



That is to say, the use of Minkowski's stress tensor is 
expected not to yield the correct force in general, whereas 
the use of T, which is formally the same as the stress 
tensor in microscopic electrodynamics, is always correct. 
Let 



P = Pind " 
M = M ind 



Pn, 

-M n 



(20) 
(21) 



be the decompositions of the polarization and the mag- 
netization into induced parts Pj, lc i,Mi nc i and noise parts 
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Pn,Mn, where the noise parts are closely related to dis- 
sipation. Substituting in Eq. I|18|) for P and M the ex- 
pressions H2()(l and (|21|1 . respectively, we see that force 
calculations that are based on Minkowski's stress tensor 
are expected to be incorrect with respect to both the in- 
duced parts and the noise parts of the polarization and 
the magnetization in general. Clearly, if — and only if — 
the aim is to calculate the force acting on bodies that 
are placed in a free-space region, then both T and T^' 
lead to the same result. 

The idea to regard [according to Eqs. (|1(J|) — (|1 3f> ] the 
Lorentz force acting on the totality of charges and cur- 
rents belonging to a system under consideration as the 
fundamental quantity is neither new [111 l2ll l22l |23j nor 
particularly hard to agree with. Despite this, the use 
of Minkowski's stress tensor or related quantities has 
still been common in the calculation of electromagnetic 
forces. In this context, let us make a few general remarks. 
The momentum that may be introduced on the basis of 
Eq. is related to the Noether symmetry expressing 
homogeneity of space. It must be distinguished from the 
pseudomomentum related to (strict) homogeneity of the 
material. In connection with the so-called Minkowski- 
Abraham controversy, Refs. 0, 0] analyze in a La- 
grangian framework the meaning of different momentum- 
like quantities by consideration of explicit (classical) dy- 
namical models of a homogeneous dielectric. In Ref. |21| . 
the homogeneous dielectric is assumed to be lossless and 
treated in some multipolar, long- wavelength approxima- 
tion (for an inclusion of magnetic properties, see pljl. 
In Ref. (24|, the homogeneous dielectric is described by 
a single-resonance Drude-Lorentz model. All the calcu- 
lations show that Eq. [together with Eqs. (|10|l and 
(|12|l ] is really the momentum balance of the macroscopic 
electromagnetic field. 



III. CASIMIR FORCE ON BODIES EMBEDDED 
IN DISPERSING AND ABSORBING MEDIA 

In classical electrodynamics, electrically neutral mate- 
rial bodies at zero temperature which do not carry a per- 
manent polarization and/or magnetization are not sub- 
ject to a Lorentz force in the absence of external elec- 
tromagnetic fields. As known, the situation changes in 
quantum electrodynamics, since the vacuum fluctuations 
of the electromagnetic field can give rise to a nonvanish- 
ing Lorentz force — the Casimir force. Its experimental 
demonstration has therefore been regarded as a confir- 
mation of quantum theory. 

To translate the classical formulas given in Sec. ITT1 into 
the language of quantum theory, let us consider linear, 
inhomogeneous media that locally respond to the elec- 
tromagnetic field and can thus be characterized by a 
spatially varying complex permittivity e(r, w) and a spa- 
tially varying complex permeability /i(r,w). Following 
Ref. [26j , we may write the medium-assisted electric and 



induction field operators in the form of 

/>oo 

E(r)=/ dwE(r,w)+H.c, (22) 
Jo 

/>00 

B(r) = / <k)B(r,w)+H.c, (23) 
Jo 

where 

E(r, lj) = iiiou J dV G(r, r', w)] N (r', to), (24) 
B(r, w) = A^oV x J dV G(r, r', c)] N (r', w). (25) 

Here, G(r, r',uj) is the classical Green tensor, which has 
to be determined from the equation 

u) 2 

V x k(t, ojW x G(r, r', uj) T e(r, w)G(r, r',u) = S(r-r') 

c A 

together with the boundary condition at infinity, and 
j N ( r ; w ) is defined by 

j N (r,u) = -wP N (r,u)+VxM N (r,4 (27) 

where P N (r, uj) and M N (r, uj) are, respectively, the (fluc- 
tuating) noise parts of the polarization P(r, uj) and the 
magnetization M(r, uj) in the frequency domain, 

P(r,w) = e [e(r, W )-l]E(r,^)+P N (r,w), (28) 
M(r,w) = Ko[l-K(r,w)]6(r,w)+^(r,w), (29) 

[n =fiQ 1 , k(t, uj) = /i _1 (r, uj)}. 

The Green tensor [as well as e(r, uj) and re(r, uj)] is holo- 
morphic in the upper uj half-plane and has the "reality" 
property 

G(r,r',-u;*) = G*(r,r',uj). (30) 

Moreover, it obeys the reciprocity relation 

G(r,r»= G T (r',r,Lu) (31) 

(the superscript T denotes matrix transposition) and the 
integral relation 

J d 3 s{[G(r,s,w)xv s ]lm K (s,w)[V s xG*(s,r',w)] 

uj 2 i 
H — 2 ^( r ) s j w)Ime(s, uj) G*(s, r', ui) > =Im G(r, r', uj), 

(32) 

where the notation 

G(r,r» x V' = -[V x G T (r,r',w)] T (33) 
has been used. 
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According to Ref. |26|. P N (r,w) and M N (r, uj) can be 
related to bosonic fields f e (r,cj) and f TO (r,o;), respec- 
tively, in such a way that the correct (equal-time) com- 
mutation relations of the electromagnetic field operators 
are satisfied, 

1 11 A 

P N (r,w) = i [fiE Ime(r, w)/7r] f e (r,w), (34) 

1 /2 A 

MnI 1 *) w ) = [-fi«oImK(r, w)/tt] 7 f m (r,w), (35) 

[/^(r.wj.^fr.w)] = 6 kl 6xx>6(r - r')6(u - u/) (36) 

(A = e,m). Note that the fx(r, w) and fj(r, cj) play the 
role of the dynamical (canonical) variables of the theory. 
Combining Eqs. with Eqs. CH, and 

yields the electric and induction fields in terms of the 
dynamical variables. 

The charge and current densities that are subject to 
the Lorentz force are given by 



J(r) 



duj p(r, w) + H.c, 
dwj(r,w) +H.c, 



(37) 
(38) 



where, according to Eqs. ©-© (p e xt = 0, j ex t = 0) to- 
gether with Eqs. l(T7|) . and 



p(r,v) = -e V{[e(r,^) - l]E(r,w)} + (iwJ^VjJr.t 



and 



(39) 



j(r,w) = -iwe [£(r,w) - l]E(r,w) 

+ VX {Ko[l-K(r,«)]6(r,«)}+j N (r l w). (40) 

Using the definitions of E(r, w), B(r, uj), and j N (r, to- 
gether with the bosonic commutation relations for the 
fundamental fields f\(r,u>) and fUr,u>), one can prove 
(Appendix El that 



0. 


(41) 


0. 


(42) 


0. 


(43) 



and 



[i fc (r),^V)] =^(r)«^(r-r'), (44) 



where the position-dependent plasma frequency f2 e (r) 
is defined by the asymptotic behavior of the permit- 
tivity for large to in the upper half-plane according to 
e(r, us) ~ 1 — r2^(r)/w 2 . The commutation relations l|41|l - 
(|44|l clearly show that p(r) and j(r) really represent mat- 
ter quantities. It is worth noting that Eq. I|44(l exactly 
corresponds to the equation obtained when — on the ba- 
sis of a microscopic description — the current density is 



explicitly specified in terms of particle velocities (Ap- 
pendix . 

If the field-matter system is in a number state [de- 
fined with respect to the number (density) operators 
fj[(r, ui){\(r, us)} such as the ground state, or an incoher- 
ent mixture of them such as a thermal state, then all 
one-time averages are evidently time- independent. Re- 
calling the bosonic character of the fundamental fields 
fx(r,w) [and fj(r,u;)] and assuming them to be excited 
in thermal states, we easily obtain, in close analogy to 
Ref. M, 



(f A (r, W )®ft,(r>')) 



coth 



2k B T 



Sxx'S(cJ-uj')S(r-r'), (45) 



(ft(r,w)®f v (rV)) 



coth 



2k B T 



1 



5xx'8(u - u')6(r - r'), (46) 



(f A (r,w)®f.y(rV)> = 0. 



(47) 



Making use of Eq. JSJJ together with Eqs. J22I) and 

we find that the correlation functions H45fl - I|47|) imply the 

correlation functions 



(L(r, W )®£(r',c')) 



j n \ i i 

h 



,2 



8{oj-u') 



coth 



2k B T 



< — \/lme(r, uj)d(r — r')^lme(r',uj') 
V x UniK^wJ^r - r')V ImK ( r ', w ')l xV'|, (48) 



(j;(r, W )®] N (r',o;')) 

h 



,2 



iN 

■6(u-u') 



coth 



2k B T 



-1 



r ^2 

< — \/lme(r, w)£(r — r')-\/lm£(r', a/) 

VWM<i(r - r')\/lmK(r',u/)] xV'|, (49) 



Vx 



and 



<L(r, W )®L(r>')> = 0. 



(50) 



Using Eqs. (|22I> and (f^SJ together with Eqs. (J23JI and 
and employing Eqs. lO, lEU), and gBl)-(G>nj>, we can 

calculate the thermal-equilibrium correlation functions of 
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the electric field and the induction field to obtain 



(E(r)®E(r')> 

= ^° J^du coth(^) Im G(r, r', a,), (51) 

(B(r) ® B(r')) 

hfX ° f °°du; cothf J^r ) V x hi) Glv. ,:•'..') • V. 



where 



(52) 



Taking the limit T — > (i.e., replacing the hyperbolic 
cotangent with unity) yields the respective ground-state 
correlation functions. Note that the correlation functions 
()51J) and l|52() inherit the reciprocity property according 
to Eq. l|3T|) . 

Now we calculate the expectation value of the Lorentz 
force [which is Hermitean — recall Eqs. (|41|l and 1|43|) ]. 



/ d 3 r( j 5E+j x B), 
Jv 



(53) 



where E and B, respectively, are defined by Eqs. I|22|) and 
(|23[l together with Eqs. (|24l) and l|25|l . and p and j, respec- 
tively, are defined by Eqs. (|37[1 and (|38[1 together with 
Eqs. i|39|) and H4U|) . Following the line suggested by clas- 
sical electrodynamics, paying proper attention to opera- 
tor symmetrization as well as regularization, and taking 
into account that the time derivative in the (quantum- 
mechanical version of) Eq. (|14|l does not contribute to 
the force, we find ( Appendix lB)l that Eqs. ifTBl and l(TB|) 
apply, where the (time-independent) stress tensor can be 
obtained, in agreement with the classical Eq. Q12[l . from 
the quantum-mechanical expectation value 

T(r,r') 

= e (E(r) ® E(r')> + ^ X <B(r) ® B(r')) 
-\ 1 [e (E(r)E(r')) + ^ '(B(r)B(r'))] (54) 

in the limit r' — > r, where divergent bulk contributions 
arc to be removed before taking the limit [recall the re- 
mark below Eq. (|16|l ] . This is always possible if the body 
under study is embedded in a material environment that 
is homogeneous at least in the vicinity of the body. If 
this is not the case, special care and additional consider- 
ations are necessary, and it may happen that physically 
interpretable results can hardly be extracted. Note that 
in the calculation of the surface integral in Eq. H15|l the 
"outer" values of the integrand should be used if dV is 
the interface between an inhomogeneous body embedded 
in a homogeneous environment (see Appendix IB|l . 

Inserting Eqs. ijSTJl and into Eq. finally yields 
the stress tensor as 

T(r, r) = lim [0(r, r') - \l Tr 0(r, r')l , (55) 

r' — >r 



■ Im G(r, r', w) - V x Im G(r, r', u>) x V' 



(56) 



As expected, the permittivity e(r, u>) and the permeabil- 
ity fi(r,u>) do not appear explicitly in Eq. (|56[) . but only 
via the Green tensor G(r, t',oj). Having removed diver- 
gent bulk contributions, we may take the imaginary part 
of the whole integral instead of the integrand in Eq. I|5t)|) 
and rotate the integration contour in the usual way to- 
ward the imaginary frequency axis, on which the Green 
tensor is real [recall Eq. I|30|) ]. In the zero-temperature 
limit, the result is simply 



0(r,r>) = -- \ 

T JO 



V G(r, r', if) + V x G(r, r', if) x V 



(57) 



For nonzero temperatures, a sum over the poles of the 
hyperbolic cotangent (corresponding to the Matsubara 
frequencies) arises instead. It should be mentioned that 
the zero-frequency contribution to the resulting series can 
be problematic if the expression in the square bracket in 
Eq. I|57l) has a singularity there, which is the case when 
permittivities of Drude type (exhibiting a pole at zero 
frequency) are employed. In fact, this unpleasant fea- 
ture expresses the conceptual limitations of a spatially 
local description of the material response to the electro- 
magnetic field, which disregards spatial dispersion. It is 
known that, for materials with (almost) freely movable 
charge carriers, this can become an issue especially at 
low frequencies (large free path lengths). Extension of 
the quantization scheme to nonlocally responding mate- 
rials would render it possible to include such materials in 
the calculation of Casimir forces in a consistent way. 



IV. HARMONIC-OSCILLATOR MEDIUM 

It is maybe illustrative to make contact with micro- 
scopic approaches to the problem. The simplest and 
most widely used model for describing linearly polar- 
izable media is quite certainly the harmonic-oscillator 
model (inclusion of magnetic properties into the model 
is still scarce). To account for dissipation, the medium 
oscillators that are relevant to the linear interaction with 
the electromagnetic field — shortly referred to as medium 
oscillators — are also linearly coupled to (infinitely many) 
heat bath oscillators (e.g., phonon modes). The effect 
of the heat bath can then be adequately taken into ac- 
count by including friction terms and associated noise 
forces in the equations of motion of the medium oscilla- 
tors. On a coarse-grained time scale, the friction terms 
are commonly regarded as being local in time (Markov 
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approximation), so that they can be characterized by 
simple damping constants. It should be noted that the 
requirement of limited time resolution implies that differ- 
ent, not strictly equivalent noise forces are acceptable in 
that regime (for details of dampin g th eory and oscillator 
models, see, e.g., Refs. [13, 113,^3 ) - 

In the context of the one-dimensional theory of 
the Casimir force on absorbing bodies, the harmonic- 
oscillator model has been used to study the interaction 
of damped medium oscillators with the transverse part 
of the quantized one-dimensional electromagnetic field, 
with special emphasis on homogeneous media |32|. Ex- 
tending the one-dimensional theory to three dimensions, 
we begin with the Heisenberg equations of motion of the 
system in the form of 

p(r, t) = -mwps(r, t) ~ m~/s(r, t) + eE(r, t) + F N (r, i), 

(58) 



S(r,t) =p(r, t)/m, (59) 



V x B(r,t)- i =Moj(r,*), (60) 



VxE(r,t) = -B(r,t), (61) 

where s(r,t) and p(r, t) are, respectively, the coordinate 
field and the momentum field of the medium oscillators, 
and 



j(r,t) = ei7(r)s(r,t) 



(62) 



is the (model) current [f?(r), number density of the me- 
dium oscillators]. Further, Fn(i", t) is the Langevin 
noise force acting on the damped harmonic oscillators 
(7, damping constant). In addition to the equations of 
motion Eqs. I|58|) (|6H and the definition (|62f> . the com- 
mutation relations J28j 



[s(r,i),p(r',i)] = — a(r-r') 



and 

[F N (r,i),F N (r',0 



2 7717 
r)(r) 



r](r) 



ihS(r 



(63) 



together with the standard commutators of the electro- 
magnetic field are required to specify the model. The 
commutator (I64f> ensures that Eq. I|t)3l) is preserved in 
time. Note that the proof given in Ref. 28] extends to the 
inhomogeneous case rj = rj(r) considered here. It should 
be stressed that the number density rj(r) is not allowed to 
have zeros (nor infinities), otherwise Eqs. (|63|) and i|tj4|) 
were not well-defined. Further, note that 



is the full electric field consisting of both longitudinal and 
transverse parts [38| , 



E«(- L )(r > *)= / dV<* ll(±) (r-r')E(r',<). 



(66) 



The transverse part may be associated with a vector po- 
tential in the Coulomb gauge and expanded into orthogo- 
nal modes in the usual way. By contrast, the longitudinal 
part is not a dynamical electromagnetic field variable but 
is (nonlocally) determined by the oscillator field as 

EM(r,t) = -e[»j(r)S(r,t)]ll/eo, (67) 

implying the conserved (model) charge density 

p(r,t) = -eV[/)(r)s(r,t)], (68) 

which is consistent with Eqs. I|60|) and (|62|l . Needless 
to say, Eqs. (|62|) and I|tj8fl do not actually represent 
the sources on a truly microscopic level but rather on 
a mesoscopic one, since the term "continuously varying 
field" applied to matter consisting of well-distinguishable 
constituents already indicates some averaging. However, 
complying with established terminology, we refer to this 
mesoscopic description as being microscopic throughout 
the paper. Note that ab initio calculations on a truly 
microscopic level would lead to time-ordered products 
in the treatment of the complicated interaction problem 
even in linear electrodynamics, because of the interaction 
with the dissipative system. However, if this interaction 
is treated in Born and (quasi-)Markov approximations, 
then the closed equations derived in this way no longer 
contain any time-ordered products. 

As the system evolves toward its dressed ground state 
as t — > 00, the model can be shown (Appendix [UJ to 
lead to the equal-time electromagnetic field correlation 
functions 

lim (E(r, t) <g> E(r'. t)) = ^ f°°6u u? Im G*(r, r', u) 

(69) 

and 

lim (B(r,t)®B(r',t)) 

t—>oo 

= -0^1 [°°Au V x Im G(r, r', w) X V' (70) 

7T JO 

if the heat bath that interacts with the medium os- 
cillators is assumed to have zero temperature. Here, 
G(r,r',u>) is the Green tensor that is the solution to 
Eq. |J2EJl, with n(r,u>) and e(r,w), respectively, being the 
model-specific quantities K(r, u) = 1 and 



e(r,uj) = 1 + 



e 2 ij(r) 



1 



Eom uJq — uj 2 — ijuj 



(71) 



E(r,i) =E"(r,i)+E- L (r,t) 



(65) 



Obviously, Eqs. I|t)9|) and l|70|) . which directly follow from 
the microscopic model under consideration, correspond 
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exactly to Eqs. (|51|l and 1(52(1 in the zero-temperature 
limit. 

The (steady-state) Lorentz force acting on the 
harmonic-oscillator matter in some space region V is 
given by 

F= lim f d 3 r(p(r,t)E(r,t)+j(r,t)xB(r,t)), (72) 

with p(r,t) and j(r, t) from Eqs. (|6"5j) and Ip^jl. respec- 
tively. At this stage it not difficult to see that the pro- 
cedure outlined in Appendix [B] yields F in the form of 
Eq. I |l 5(1 together with Eq. 1(16(1 . where the stress tensor 
has exactly the form of Eq. 1(55(1 together with Eq. ((57(1 . 
This result shows that the microscopic approach to the 
Casimir force fully confirms the macroscopic approach 
as given in Sec. IIIII where the calculations were based 
on the quantized macroscopic electromagnetic field, with 
the matter phenomenologically described in terms of 
Kramcrs-Kronig consistent response functions. Thus, the 
Casimir force acting on a macroscopic piece of matter 
may be viewed as "just" the (quantum) Lorentz force on 
the constituting charges and currents, which, in a macro- 
scopic description, can be expressed in terms of the (in- 
duced and noise) polarization and magnetization — a con- 
ceptually straightforward and satisfactory point of view. 



V. CASIMIR FORCE IN PLANAR 
STRUCTURES 

Let us apply the theory to a planar magnetodielectric 
structure defined according to 



e(r,uj) = 



e-(z,uj), 

Ej(w), 

e + {z,uj), 



z < 0, 

< z < dj, 

z > dj, 



z < 0, 
< z < dj 
z > dj. 



(73) 



(74) 



To determine the Casimir stress in the interspace 
< z < dj , we need the Green tensor in Eq. 15611 for both 
spatial arguments within the interspace (0 < z = z' < dj ) . 
The Green tensor is well known and can be taken, e.g., 
from Ref. |3^|. Since the transverse projection q of the 
wave vector is conserved and the polarizations a = s,p 
decouple, the scattering part of the Green tensor within 
the interspace can be expressed in terms of reflection co- 
efficients rj ± =r?± (u>, q) referring to reflection of waves at 
the right (+) and left (— ) wall, respectively, as seen from 
the interspace. Explicit (recurrence) expressions for the 
reflection coefficients are available if the walls are mul- 
tislab magnetodielectrics like Bragg mirrors [271 l33l l34| . 
(For continuous wall profiles, Riccati-type equations have 
to be solved [33. ) In the simplest case of two homo- 
geneous, semi-infinite walls, the coefficients rj ± reduce 



to the well-known Fresnel amplitudes. In the case first 
treated by Lifshitz [35| . the interspace is empty and the 
walls are nonmagnetic. 



A. Casimir stress within a nonempty interspace 

For the sake of generality, we first leave the wall struc- 
ture unspecified. By modifying the expression for the 
scattering part of the Green tensor given in Ref. [34[ to 
also account for magnetic properties, from Eq. 1(55(1 to- 
gether with Eq. (|5rJ)l (without the bulk part of the Green 
tensor) it then follows that the relevant stress tensor el- 
ement T zz (r, r) in the interspace < z < d j can be given 
in the form of 



fj,j(0j) 



x Re / dqq 



Pj(u,Q) 



9j(z,u3,q) (75) 



(q= |q|), where the function gj(z,u>,q), which in general 
depends on the position z within the interspace, reads 

9j(z,iv,q) 

= 2[#(1 + nj 2 ) -q 2 (l- nf)]Djy j+ r S j_e^ 

+ 2[#(1 + nf) + q 2 (l n^)]Djy j+ r^e 2 ^ 

- (13 2 + q 2 )(l - nJ 2 )D]f [r]_e 2l ^ z + r s j+ e 2i0 ^-^] 

+ (p 2 + q 2 ){\ - nJ 2 )Dji y_e 2 ^ + fP + e 2 ^(*-»)] , 

(76) 



with the definitions 

n 2 = n 2 (uj) = ej{u)nj(u), 
/3^/3 J (c,g) = (^V 2 /c 2 - 9 2 ) 1/2 , 



D j<r = D ja (cj,q) = 1 



(77) 
(78) 
(79) 



Note that the equations Dj a (uj,q) = determine, for 
real q, the frequencies of the guided waves in the pla- 
nar structure, which are of major interest in all "mode 
summation" approaches. (In the presence of material 
losses, however, these waves have complex frequencies 
and are not ordinary normal modes.) For practical rea- 
sons, it may be advantageous to transform the integral 
over real frequencies in Eq. ((75(1 into an integral along 
the imaginary frequency axis by means of contour inte- 
gral techniques [cf. Eqs. I(55|) and ((57(1 ], In particular, in 
the zero-temperature limit, Eq. ((75(1 may be rewritten as 



T„(r,r) = — j / d£ 
8tH Jo 



From the derivation it is obvious that the stress for- 
mula ((75(1 [together with Eq. 1(76(1 ] allows for an arbi- 
trary linear, causal interspace medium. By contrast, 
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Minkowski's stress tensor 
fa = 1) 



[Eq. 03] leads to H 



T (M) (r)] 



h 

'2^2 



, ( fit" 



Re dqqp.Y" 



'j+'j- 



(81) 



From Eq. Q76[l it is easily seen that for an empty inter- 
space, i.e., Ej = jij = 1, gj(z,uj,q) becomes independent 



of z and simplifies to 



gj(z,u,q) -> 9j(u},q) 



4$ 



T — S,p 



D, 



(82) 



In this case, and only in this case, Eq. i[75J) reduces to 
Eq. (|81(l , from which in the case of semi- infinite (homoge- 
neous) dielectric walls Lifshitz's well-known formula |35j 
can be recovered. As already mentioned, formulas of the 
type of Eq. I|81|) [which need not necessarily be derived 
within the stress tensor approach to the Casimir force] 
have been claimed to apply also to the case where the in- 
terspace is filled with dielectric material 0,0, at least if 
the material is nonabsorbing Q (see also the textbooks 

00 and references therein). Since (r,r) docs 
not depend on the position z within the interspace, ap- 
plication of Eq. I|81(l implies the very paradoxical result 
that the force acting on any slice of material selected 
within the interspace vanishes identically, regardless of 
the presence and arrangement of the remaining material 
(in particular, regardless of the yet unspecified walls). 
This unphysical result clearly shows that Eq. I|81|l can- 
not be valid if the interspace is not empty, not even if it 
may be justified to regard the interspace medium as non- 
absorbing. In contrast, the stress T zz (r, r) obtained from 
Eq. I|75l) [together with Eq. J7SJ)] is not uniform within 
an interspace if the interspace is filled with a medium. 
Hence it gives rise, in general, to a nonvanishing force on 
a slice of interspace material, and no paradox appears. 

Let us return to the stress formula (|75|l [together with 
Eq. iJZB}]. It is not difficult to see that, for a nonempty 
interspace, the q integral in Eq. I|75|l fails to converge at 
z = and z = dj, i.e., on the interfaces where the differ- 
ent materials are in immediate contact with each other. 
Mathematically, the reason for this divergence can be 
seen in the fact that the reflection coefficients obtained 
under the assumption of infinite lateral extension of the 
system do not approach zero as q tends to infinity. How- 
ever, large values of q correspond to very oblique trav- 
eling waves. In any real planar setup of finite lateral 
extension, such high-g waves clearly do not contribute to 
the q integral at all; they are not reflected but walk off 
instead. Note that a divergence of exactly the same type 
already appears also in the standard case of an empty in- 
terspace in the limit dj — » 0. In order to (approximately) 
take into account the finite lateral extension of an ac- 
tual planar setup, an appropriately chosen cutoff value 



(depending on the lateral system size) for the reflection 
coefficients at high q values could be introduced, thereby 
rendering the q integral finite. Of course, a more satis- 
factory approach would be to abandon the translational 
invariance from the outset, which, however, leads to se- 
rious mathematical difficulties since waves with different 
polarizations and transverse wave vectors are then no 
longer decoupled. 

Since, according to Eq. i(T5|) . the Casimir force acting 
on a body is given by the integral of the stress tensor 
over the surface enclosing the body, the stress tensor on 
its own is of less importance. What is really important is 
the integral force value over a closed surface. To obtain 
the force (per unit area) acting on a (multilayered) plate 
of infinite lateral extension, the stress on the two sides of 
the plate must be taken into account. As the example in 
Sec. IV Bl shows. it may then happen that the parts of the 
stress tensor that diverge when the plate is approached 
from the two sides cancel each other out. In such a case, 
the Casimir force (per unit area) on a plate remains well 
defined even if its lateral extension is assumed to be in- 
finite. 



B. Casimir force on a plate in a nonempty cavity 

In order to make contact with recent work on the 
Casimir force on bodies embedded in media [4j, let us 
calculate the force acting at zero temperature on a ho- 
mogeneous plate in a nonempty planar cavity, according 

z-direction 








FIG. 1: Homogeneous plate embedded in a nonempty cavity. 
The cavity medium on the right and left sides of the plate is 
the same. 

to the five-region setup as sketched in Fig. ^ The cavity 
walls are labeled by I — and I = 4, the plate by 1 = 2, 
and the cavity regions that are filled with the medium 
the plate is embedded in are labeled by I — 1 and / = 3, 
with e(oj) = ei{bj)=e^{ui) and /u(w) = /zi(u;)=/Z3(w). The 
total (volume) force per unit transverse area acting on 
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the plate can be obtained by (vectorial) addition of the 
two force contributions from the two sides of the plate. 
Application of Eq. 1(80(1 then yields the total force per 
unit transverse area in the form of 



-gi(di,i£,q)] (83) 

Mi (ie,g) = /%(i£,g)]. 
For a quantitative comparison with specific results ob- 
tained in Ref. 4] on the basis of Minkowski's stress 
tensor, we make the following simplifying assumptions. 
We assume that (i) all the reflection coefficients can be 
regarded as being almost constant, and (ii) the reflec- 
tion coefficients rf_|_ and r 3 _ can be approximated by 
the (same) single-interface (Fresnel) reflection coefficient 

1/2- 



From an inspection of Eq. 1(86(1 it is seen that there is 
no divergence; the integrals are well behaved. It is worth 
noting that even without application of the approxima- 
tion scheme, the integrals in the basic formula 1(83(1 do not 
diverge. The reason is that, for a chosen value of f , the 
coefficients r 3 _ (if, q) and rf + (zf , q) tend exponentially to 
the same single-interface Fresnel coefficient r^, 2 (i^,q) as 
q goes to infinity, as may be seen from relations like 



1/2 



+ e 



2i0 2 d 2 



'2+ 



' 1+ 



'l/2 e '2+ 



'3/2 



= 2i/3 2 <i 2 



'3/2 C '2- 



1/2 



'3/2 



if q 



if q 



(87) 



(88) 



Physically, these assumptions mean that (i) the together with the relation 



3/2 '1/2 



(valid for arbitrary 



distances d\ and ^3 between the plate and the cavity 
walls must not be too small, and (ii) the plate must be 
thick enough. Moreover, the approximation scheme im- 
plies that the permittivity and the permeability of the 
medium the plate is embedded in can be replaced with 
their static values briefly referred to as e and [i in the 
following, with n = a/e/Z being the static refractive in- 
dex. From Eq. 1(76(1 it then follows that the difference 
of the functions 173(0, if, q) and g\ (d\, if, q) appearing in 
Eq. i(83|) can be approximated according to 



#3(0, if,?) - 9i{di,i£,q) 
1 1 



(T—S,p 



D. 



3ct 



D 



la 



A a q 2 1-- 

n z 



1/2 



A.(/3 2 



1/2 



1-- 

n 



r a p 2i0d 1 



(A CT = 5 ap - Sera), where 



tr 2if3d 3 er 2ifld\ 



r\_e 



1 - D 3a 1 - D l0 



(84) 



D 



'An 



(85) 



Substituting Eq. (JSIJ together with Eq. (JHSJl into 
Eq. H83|) . we (approximately) obtain 



values of f and q). Note that ifo — > —00 if q — > 00. As 
a consequence, the divergent contributions to the q inte- 
gral in Eq. 183(1. which would arise from ^3(0, if, q) and 
gi(di,i£ l , q) separately, combine in a convergent fashion. 
Thus, for the setup under study, a q cutoff need not be 
introduced. 

Let us return to Eq. 1)86(1 . If the two walls and the 



plate are almost perfectly reflecting, i.e. 



'3+- 



1/2- 



'. A CT , then standard evaluation of the integrals leads 



tO (n = y / £/l) 



F = 



HCTT 2 

240 



3 + 3e^J \dj d\ 



(89) 



In particular, if only one wall is present, say the left one, 
then Eq. ((89(1 reduces to {d 3 ^00, d\=d) 



F 



HcTT 

240 



1 

3e/i 



1 



(90) 



which is the generalization of Casimir's well known for- 
mula [33 for the force between two almost perfectly re- 
flecting plates separated by vacuum [fi=e= 1 in Eq. 190(1 ] 
to the case where the interspace between the plates is 
filled with a medium of static permeability \x and static 
permittivity e. 

In order to compare Eq. 1(89|) with the force formula 
obtained on the basis of Minkowski's stress tensor, we 
note that the use of Minkowski's stress tensor for a non- 
magnetic medium leads to [see Eqs. (3.6) and (3.7) in 
Ref. 0] 



87r Jo Jo iP 



2/3 2 (l + 4l -AX(n 2 - iWr 



1 1 

1 



1/2 



T/2 



2A CT <? 2 1 - 



1 



-2i/3di _ 1 



(91) 



(86) in place of Eq. 1(86(1 with /1 = 1. For an almost perfectly 
reflecting plate in a cavity with almost perfectly reflecting 
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FIG. 2: (a) The Casimir force F given by Eq. fg£ (solid 
curve) is shown as a function of the medium refractive index 
n — ^/e (ji=i) for chosen distances di and d%. For comparison, 
7^( M ) given by Eq. 1921 (dashed curve) is shown, (b) The ratio 
p( M ) J p j s shown as a function of the medium refractive index. 



walls, standard evaluation of the integrals in Eq. 191(1 
then yields, in place of Eq. (|S9^I . 



pw 



HcTT 1 



240 \4 df, 
which in the limit d$ — > 00 reduces to (di = d) 

p(u) = J_ J_ 

240 y/e d 4 ' 



(92) 



(93) 



Note that Eq. I)98|l corresponds to the result derived in 
Ref. [33] by means of mode summation methods. Com- 
paring Eq. (HO) with Eq. © [or Eq. ^ with Eq. J33J)], 
we see that 



\p\ < |^(M)| ; 



(94) 



i.e., the absolute value of the force is (n > 1) always 
smaller than that predicted from Minkowski's stress ten- 
sor. Introduction of a (polarizable) medium into the 
interspace is obviously associated with some screening 
of the plate, thereby reducing the force acting on it. 
Since the internal charges and currents of the interspace 



medium are not fully taken into account in a theory that 
is based on Minkowski's stress tensor or an equivalent 
formalism, the screening effect is underestimated and 
consequently the force calculated in this way is overes- 
timated. Although the assumptions made to derive the 
results given above are rather restrictive, the comparison 
of Eq. <|%9"|) with Eq. (|9"2"|) clearly shows that the correct 
inclusion of the medium into the theory can give rise to 
noticeable effects (see Fig- EJ ■ 

A consequence of the approximation scheme employed 
in this section is the appearance of the real values of the 
static permittivity and the static permeability of the in- 
terspace material in Eq. 189|) . However, the basic equa- 
tion (|83l) is of course valid for arbitrary linear magne- 
todiclectric media with Kramers-Kronig consistent per- 
mittivities and permeabilities. The influence of material 
dispersion and absorption comes into play when the dis- 
tances di and/or d^ are decreased. The behavior of the 
permeability and the permittivity at nonzero frequencies 
becomes then important. 



VI. SUMMARY AND CONCLUSIONS 

On the basis of (i) the quantized macroscopic electro- 
magnetic field in the presence of causal linear magne- 
todielectric media without spatial dispersion and (ii) the 
Lorentz force acting on the internal charges and currents 
of the medium, we have derived general expressions for 
the Casimir force acting on magnetodielectric bodies em- 
bedded in a common magnetodielectric medium. All the 
matter has been allowed for being dispersing and absorb- 
ing. Specializing to planar structures, we have general- 
ized Lifshitz-type formulas (being valid for empty inter- 
spaces) to the case where the interspaces are filled with a 
magnetodielectric medium. In this context, we have ana- 
lyzed the failure implied by basing the calculation of the 
Casimir force on Minkowski's stress tensor — a method 
that has been used in the literature but has never been 
proven correct. Interestingly, Lifshitz himself did not ad- 
dress nonempty interspaces in his seminal article |35|. 

For comparison reasons, we have studied in some detail 
the Casimir force acting on a homogeneous plate embed- 
ded in a medium in a planar cavity. Applying standard 
approximations such as high reflection, we have explic- 
itly demonstrated that when the plate is embedded in 
a medium, then the force can noticeably differ from the 
result obtained on the basis of Minkowski's stress ten- 
sor. By the way, we have given the correct extension of 
Casimir's original formula for the force between two per- 
fectly reflecting plates to the case where the interspace 
between the plates is filled with a medium. 

In order to make contact with microscopic theories, we 
have also described the matter microscopically, by em- 
ploying the model of damped harmonic oscillators, which 
is widely used for treating dielectric matter. Solving the 
quantum-mechanical equations of motion of the overall 
system (with the heat bath assumed in its ground state), 
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we have calculated the Lorentz force acting on a cho- 
sen matter element. The result obtained in this way ex- 
actly corresponds to the general result obtained from the 
macroscopic approach. This clearly shows that the use of 
Minkowski's stress tensor to calculate the Casimir force 
is wrong in general, even if the matter may be regarded 
as being nonabsorbing. 

Note added. Instead of Eq. I|84|l . it may be advantageous 
to use the exact equation 



53(0, i£,g) - gi{d!,i£,q) 



£ 2 



P 2 1 + 



A ct9 2 1 



o2 ,a2 



t aZ ) 



r 3+ e 



N° 



(95) 



where 



N a = 1 - r CT (r^e 2i/3dl + r% + e 2i(3d3 ) 

+ {r a2 - t a2 )r1_rZ + e 2i ^ dl+d3 \ (96) 

with r a = rf/g = r% , x and t a = t^, 3 — ig /j being single-plate 
reflection and transmission coefficients, respectively. We 
thank Marin-Slobodan Tomas for this suggestion. 
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APPENDIX A: PROOF OF EQS. j4lT>-(f44l 

Using Eqs. Eqs. J23J), P5)l- and (J2BJI, we express p(r, us) 
and j(r, u) as defined by Eqs. (f?TT?f) and respectively, 
in terms of J N (r, oj) to obtain 

p(r, w) = ^ V / dV G(r, r', c)L(r', u ), (Al) 



j(r,w) = V x V x 



A'G(r,r»j N (r', U ). 

(A2) 



By combining Eqs. (2ZJ, JS}, and (J3HJ) with the stan- 
dard bosonic commutation relations for the fundamental 



fields i\(r,u>) and f|(r, w), it is not difficult to show that 
j N (r, u>) and j + (r, lo) obey the commutation relation 

[l Nk ( r --)'l^-0] = ^%--') 



— \/lme(r,w) <$(r — r')\/lm£(r', 1 

c 2 



V x ^/lm/t(r,a;) £(r — r')\/lmK(r', a/) x V' 



(A3) 



From Eqs. ijHJ, (J2SJ, (|A1|) . and (|A2ft together with the 
commutation relation (|A3|) . we derive, on recalling the 
Green-tensor relations (|31|) and l|32|) . 



[p(r, «), E f (r', a/)] = *(w - w') V Im G(r, r', 



2 

7T C 



■[pt(r, W ),E(r',u/)], (A4) 



^(r, W ),S;(r',a;')] 



7T 



-[jt(r,a,),^(r',a/)], 
[pM.bVV)] 



VxVx 5- ImG(r,r',w)x V 



k! 



(A5) 



a/1 V Im G (r, r', w) x V 

7T C 2 
t/ 



[>t(r,w),B(r',a/)], 



(A6) 



[j fc (r,u,),^VV)] 
= (5(tj — u/) 



Vx Vx llm G^(r, r',w) 



A7 



= [|;(r,a;) ) ^(r',a;0], (A7) 
where 

G- L (r,r',cj)= /" d 3 s G(t,s,lj)8 ± {s - t'). (A8) 

Note that in Eq. i|A6(l . G ± (r, r', w) may be replaced with 
G(r,r',cj), because of the operation x V. 



Equations (|A4|) and (|A5I) obviously imply the commu- 
tation relations 

[p(r), E(r')] = jf " do; jf °° do;' j [p(r, w ), E 1 (r', «')] 

+ [pt(r, W ),E(r',u/)]}=0 (A9) 



and 



(r), Bi (r')] = jf °° dc jf °° da/ j [j_ k (r, u), B* (r', u')] 
+ ^(r,w),^(r>')]}=0, (A10) 
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and hence Eqs. (|41|l and l|43[l are seen to hold. 
Note in particular that the commutation relation 
[p(r),E ± (r')] = is valid. From Eqs. JSfjJ and ifATf . 

respectively, it follows that 

[p(r),B(r')] = j / dwwVImG x (r,r» x V' 

71-0 Jo 

(All) 

and 



[j fe (r),£^(r')] 

7T 



dw W 



V x V x 



U3 



ImG^r, r',oj) 



kl 



(A12) 



To further evaluate the integrals in Eqs. (|A11|) and 
we recall the asymptotic behavior of e(r, oj) and 
n(r,u) for large to in the upper half-plane, viz., 



e(r, u) ~ 1 
K(r, w) ~ 1 



fi|(r) 



(A13) 
(A14) 



Substituting Eqs. (|AT3)l and (|Al4|l into Eq. (jSHJ), we eas- 
ily see that the Green tensor asymptotically behaves like 



G(r,r» 



<S(r-r') 



(A15) 



for large 10 in the upper half-plane. Thus, on recalling 
Eq. (|30|) and the holomorphic behavior of the Green ten- 
sor, we may evaluate the integral in Eq. (|A1 l|l to prove 
Eq. 



[/5(r),B(r')] 



P / dwwVG 1 ^^'^) x V 



7TC 



/ dwwVG x (r,r',w) x V 



= -tW« i (r-r')xV' = (A16) 

(P denotes principal value). Here, we have replaced the 
principal value integral along the real frequency axis by a 
contour (C) integral over an infinitely large semicircle in 
the upper half-plane and have used Eq. (|A15(I . Note that 
there is no extra pole contribution from uo = (2|| . To 
evaluate Eq. I|A12I) . we take into account that, according 
to Eq. l(2*?j|) . the relation 



V x V x - — ) Im G ± (r,r / ,uj) 

,2 



Im — [e(r,^)-l]G x (r,r» 



V x [1 — K(r,w)]V x G ± (r,r' ,lo) > (A17) 



may be used on the real u> axis. Inserting this relation 
into Eq. (|A12|I and recalling general properties of e(r, oj) 
and «;(r, w), we see that the evaluation of Eq. (|A12J| can 
be done in exactly the same way as the evaluation of 
Eq. (|A11|) . Thus, making use of Eqs. l|A13jl . (|A14|I . and 

jus), we derive 



[i,(r),^(r')] 



L 1 dww^[£(r,w)-l]Gj!i(r,r» 



= ^ 2 (r)<^(r-r'), 



(A18) 



which is Eq. ifUfl. 

For a consistency check of the commutation relation 
(|A1S|1 . let us consider a set of atoms, with each of them 
having one valence electron (e, charge; m, mass). Let r A 
be the (fixed) positions and the relative coordinates of 
the electrons. The microscopic (electron) current density 
is then given by 



j( r ) = e^2s A S(r -r A - s A ). 



(A19) 



By assuming minimal coupling and Coulomb gauge, the 
canonical momenta of the electrons commute with the 
vector potential A(r), whose conjugate momentum field 
is — £oE ± (r). Hence, we derive 

[j fc (r),4V)] 

S(r -r A - s A ) [A k (r A + s A ), E^{r')] 



e 
m 



ih e 
Eq m 



^5{T-r A -s A )8^{v A +s A -r') 



in e x -> 



5(r - r A - s A )6£i(r - r') 



(A20) 



In the macroscopic theory, the sum of the 5 functions in 
Eq. i|A20|) is expected to be replaced according to 

S(r-r A -s A )^J2 A ( r - r A-s A ), (A21) 



A 



where A(r) is a well-behaved function with unit integral, 
J d r A(r) = 1. Further, in order to produce reasonable 
coarse-graining, A(r) must be sufficiently flat so that the 
change of A(r) on atomic length scales can be regarded 
as being negligibly small. With the acting on well 
localized electronic bound states, we may hence write 



A(r - r A - s A ) ~ A(r - r A ). 



(A22) 



Thus, 



^A(r-r A -sU)~^A(r-r A )=?7(r), (A23) 

A A 

where rj(r) is the number density r/(r) of the atoms, and 
the macroscopic version of Eq. <|A20|I reads 

[j k (r),EHr')]=--v(r)5£ l (r-r>). (A24) 

L J Eq m 
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From a comparison of Eq. i|A24(l with Eq. QA18jl . the 
relation 



e 2 ri(r) 



(A25) 



is suggested to be valid, which is in full agreement 
with the harmonic-oscillator model permittivity given by 
Eq. GUI- 



APPENDIX B: QUANTUM LORENTZ FORCE 

Using Maxwell's equations Q~0 (promoted to oper- 
ator equations) together with the commutation relations 
(|41|l and l|43l) and relations of the type 

E(r)x[V x E(r)] =-[Vx E(r)] x E(r) 



±E(r)E(r) - E(r) ® E(r) 
- [VE(r)]E(r), 



(Bl) 



we derive 



p(r)E(r)+j(r) x B(r) 

- [VT(r,r')] r , =r -[V'T(r,r')] r , = 



|[B(r')xE(r)] r , =r , 



Of 



[E(r)xB(r')] r , =r , 



(B2) 



where 



T(r,r') = e [E(r) ® E(r') -\l E(r)E(r')] 

+ [B(r) ® B(r') - \l B(r)B(r')] (B3) 

is a reciprocal operator function of two spatial variables, 



T(r,r')= T T (r',r), 
because of the commutation relations 

[E(r),E(r')]=0=[B(r),B(r')]. 



(B4) 



(B5) 



Since the left-hand side of Eq. I|B2|) is Hcrmitian, so is ei- 
ther of the two alternative right-hand sides, which means 
that symmctrization is not necessary. Thus, Eq. (|llfl is 
also valid as an operator equation, and the steady-state 
equations l|15|) and (|16fl apply, with the stress tensor be- 
ing defined by Eq. I|54fl in the limit r' — > r. 

To perform the limit, we write the force acting on some 
space region V in the form of 



lim 



d 3 r / dV (5 e (r-r') 



|[VT(r,r')] + [V'T(r,r')]}, (B6) 



where 6 e (r — r') approaches 8(r — r') as e tends to zero. 
For instance, one could choose 



^(r) = (4^ 2 )- 1 <5(|r|-e), 



(B7) 



which corresponds to an average over a spherical surface 
of radius e. Let us first consider the case in which the 
material properties are homogeneous everywhere, except 
at the surface of the volume V , where they may change 
abruptly. The function T(r, r') can then be uniquely 
decomposed into a bulk part, which is divergent at r' = r, 
and a scattering part, which is well behaved at r' = r, and 
we may write 

V T (scat) (r,r) 

= [VT( scat )(r,r')] r , =r +[VT( scat )(r,r')] r , =r . (B8) 

For the scattering part, the limit e — * simply restores 
the 5 function, so Eq. (|B6|) becomes 



F= f rf'VVT (scat) (r,r) 
Jv 

+ lim J d 3 r J dV4(r-r'){[VT (bulk) (r,r')] 
+ [V'T (bulk) (r,r')]}. (B9) 



The second term on the right-hand side of Eq. (|B9(l . 
which arises from the bulk part, vanishes, as can be seen 
from the following argument |39|. Since the bulk part is 
a function of r — r', it follows that 

J d 3 r|dVi £ (r-r')Vr (bulk, (r,r')=yb(e), (BIO) 

where b is some vector that depends only on the pa- 
rameter e, and in this way selects, somewhat artificially, 
a particular direction in space. However, the bulk part 
corresponds to a setup where the whole space is filled 
with homogeneous and isotropic material, implying that 
such a preferred direction does not exist, and we can con- 
clude that lim e ^o b(e) = 0. To apply the divergence the- 
orem to the first term on the right-hand side of Eq. I|B9|) 
and transform the volume integral into a surface integral, 
we note that if the material properties change discontin- 
uously at the surface dV of the chosen volume V [cf. 
Eqs. (|2HJ) and igUJ)], then T (scat) (r, r) is also discontin- 
uous there. In view of the macroscopic description, it 
is clear that the material properties can be regarded as 
changing continuously across a sufficiently thin boundary 
layer. To include the net change across such a boundary 
layer, the "outer" values of the integrand should be taken 
(indicated by dV+), 



F= f daT (scat) (r,r). 



(Bll) 



In order to establish the validity of Eq. 1|B11(I for the 
more general case of varying material properties inside 
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the chosen space region (whose vicinity is again assumed 
to be homogeneous), one has to return to Eq. I|B6|) and 
decompose V including a thin boundary layer as de- 
scribed above into sufficiently small, nonintersecting cells 
V,. Summing over all cells, one can then show, by using 
similar arguments as above, that in the limit of vanish- 
ingly small cells, Vi — > 0, Eq. (|B11|) is obtained. 



APPENDIX C: PROOF OF EQS. (|69j> AND (f?0jl 

Combination of Eq. J5H| with Eq. yields the 

second-order differential equation 

ms(r, t) = — mw 2 s(r, t) — m7s(r, t) + eE(r, t) + Fn(i", t), 

(CI) 

and combination of Eq. I|6(J|) with Eqs. I|(jlfl and Hti2|) leads 
to 

1 d 2 - 

-2 ^E(r, ()|VxVx E(r, t) = -e M cMr)s(r, t). (C2) 

We are interested in the solution to Eqs. IjClfl and (|C2|) 
which is reached in the limit t— >oo, thereby being inde- 
pendent of the initial conditions. We may represent it in 
terms of Fourier integrals according to 



/(*) = 



/OO l. 
?e iwt /(t). 
-oo 2?r 



(C3) 

Note that the uj integrals should be treated as principal 
value integrals (with respect to ui = 0) if necessary. From 
Eqs. 1|C1(I and l|C2|) it follows that the Fourier transforms 
s(r,u;) and E(r, w) of s(r, t) and E(r,i), respectively, are 
determined by 



s(r, ui) = [muQ — mu 2 — imjui} 1 [F N (r, u>) + eE(r, w) 



(C4) 



and 



V x V x E(r, ui) E(r, ui) = e/i ?7(r)a; 2 s(r, ui). (C5) 

Substituting s(r, ui) from Eq. I|C4(I into Eq. (|C5|I and re- 
arranging, we obtain 

uj 2 

VxVxE(r,w)--£(r,w)E(r,w) = i/i u>J N (r, w), (C6) 
&■ 

where e(r,w), which is given by Eq. (|71|l . defines the 
permittivity of the harmonic-oscillator medium, and 



I N ( r ^) 



[e(r,c;)-l]F N (r, W ) 



(C7) 



is the current density associated with the Langevin force. 
The unique inversion of Eq. I|C6|1 is 

E(r, u) = i^ui f dV G(r, r', ^)J N (r', w), (C8) 



where G(r, r', ui) is the Green tensor that, for n(r, ui) = 1 
and e(r,ui) from Eq. I)7ip. solves Eq. I|26|l together with 
the boundary condition at infinity. 
To prove Eq. 1|69[) . we write 

lim (E(r,t) <x> E(r',i)) 

£ — »oo 

/OO /*00 
dw / dw'e- t(wW) *(E(r,w) ®E(r',o/)), 
-oo J-oc 

(C9) 

where, according to Eq. I|C8(1 [together with Eq. (|31|l ]. 

(E(r,w) (gi E(r',w')) = d3fi / d 3 s'G(r, s, cj) 

x (J N (s, w) <8> J N (s', c/)) G(s', r', u'). (CIO) 
If the heat bath is in the vacuum state, then 

rwyh 



(F N (r,u;)®F N (rV)) 



7rry(r) 



S(r-r') 



x / dcu" lu"5{lo" - ui)S{ui" + J) (Cll) 
Jo 

holds 28], and we find, on recalling Eq. (|C7|) . 
(J N (s,c,)®J N (s>')> - [e(s,w)-l][e(sV)-l] 
x jfe-s') /""dfc/VW'-wWf/'+a/). (C12) 



Combining Eqs. (|C9j) . (|C10|) . and (|C12jl . we derive 



lim (E(r, t) <g> E(r', i)) 



rwyh 



duiui 5 



[e(s,w)-l][e(s,-w)-l] 



x /d^Gtr,^) 1 " 1 "'^ ~ ; " J G(s,r',-u) 

(C13) 



77(s) 



From Eq. I|71|) it follows that the relation 

[e(s,cj) — l][e(s, — uj) — 1] e 2 Ime(s,o;) 



(C14) 



7?(s) £ ^7 W 

is valid for real w. Hence, we may rewrite Eq. (|C13j) as 
lim (E(r, t) cg> E(r', t)) = ^° f° dco ^ 



x y d 3 s G(r,s,w)Ime(s,w)G(s,r',-w), (C15) 

which by means of Eqs. H3()(l and (|32|) eventually leads to 
Eq. PI). 

To calculate 

lim (B(r,i) <g> B(r',t)) 

i" — s-OO 

/OO /"CO 
dw / dw / e- i ^ +w ''*<B(r,w)®B(r',w')> ) 
-oo J-oo 

(C16) 
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we express (B(r, us) ®B(r',u/)) in terms of (E(r,w) ® By means of Eqs. (|C10|I . (|CJ12|) . and l|C14|) [together 
E(r'V)), by using Eq. dHTJ) in the Fourier domain, with Eqs. ^ and it is now not difficult to prove 

Eq. H70(l . Note that there are no problems at u = 0. 

B(r, lu) = {iLu)- 1 ^ x E(r, u). (C17) 
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